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ABSTRACT 

We have calculated the ionization of low density gas cooling from a high 
0 

temperature. (->*T0“K')T~ The evolution during the cooling is assumed to be 
isochoric, isobaric, or a combination of these cases. The calculations are 
used to predict the column densities and ultraviolet line luminosities of 
highly ionized atoms in cooling gas. In a model for cooling of a hot Galactic 
corona, it is shown that the observed value of N(N V) can be produced in the 
cooling gas, while the predicted value of N(Si IV) falls short of the observed 
value by a factor of about 5. The same model predicts fluxes of ultraviolet 
emission lines that are a factor of 10 lower than the claimed detections of 
Feldman, Brune, and Henry. Predictions are made for ultraviolet lines in 
cooling flows in early-type galaxies and clusters of galaxies. It is shown 
that the column densities of interest vary over a fairly narrow range, while 
the emission line luminosities are simply proportional to the mass inflow 
rate . 
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I. INTRODUCTION 

In the past 15 years, it has become clear that hot gas (T > 10 R K) is a 

significant component of the volume of the universe, including the space 

within and between galaxies. In many cases, cooling of the gas is expected; 

these cases include gas in the galactic halo (Shapiro and Field 1976; 

Chevalier and Oegerle 1979), in spiral galaxies, in early-type galaxies 

(Forman, Jones, and Tucker 1985), and in clusters of galaxies (Fabian, Nulsen, 

and Canizares 1984). When the gas cools through the temperature range of 
4 5 

about 3 x io - 3 x 10 K, ions are present which can give strong ultraviolet 
lines (Savage 1984). Since the gas cools relatively rapidly through this 
temperature regime, thermal equilibrium is unlikely so that time-dependent 
ionization must be taken into account in making predictions for the 
ultraviolet lines. 

In this paper, we present calculations of the column densities of ions 
and line emissivities from cooling gas. The calculations treat the isochoric 
(constant density), isobaric (constant pressure), and intermediate cases. The 
results are compared with observations of galactic gas and predictions are 
made for cooling gas in galaxies and clusters of galaxies. 

II. CALCULATIONS 

We find the column densities of various observable ions in cooling gas. 

We assume that the cooling is steady and can be characterized by a cooling 

column density of H atoms per unit time, N , and an initial H density, n . We 

H O 

also assume that the cooling process is isobaric, isochoric, or is initially 
isobaric with a transition to isochoric. This should cover the likely range 
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of pressure evolution, unless there is strong evolution due to hydrodynamic 

flow. We assume that heat conduction can be neglected. Under these 

assumptions, the column density of an ion Z 1 of element Z produced between 

temperatures T . and T is 
min max 



N [ n ) max 

k ~ | — | / 

n 0 [ n Hj T . 

min 



.1 2L. 


f. dT 

i 

(n/n ) A 
o 


m 


where k is the Boltzmann constant, n/n is the abundance of the element 

L H 

relative to hydrogen, s is 1 for isobaric and 0 for isochoric evolution, x is 

the number of particles per H atom, X Q is the number of electrons per H atom, 

f. = n ./n n is the H density in the evolving gas, and A is the radiative 
1 ^ ^ z 

cooling function, such that A n„n is the energy loss rate per unit volume. 

H e 

The solution of equation (1), especially through A and f . , requires the 
calculation of the ionization of a number of elements in a cooling gas. A 
separate calculation for each type of pressure/density evolution is required. 
Kafatos (1973) and Shapiro and Moore (1976) have performed such calculations 
for an isochorical ly cooling gas. The calculations here are similar except we 
use a variety of pressure evolutions and use the atomic data of Raymond and 
Smith (3977), updated by Raymond and Smith (3984). The elements included are 
H, He, C, N, 0, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni, with abundances from Allen 
(3973). The results for isochoric cooling beginning at a temperature of 30 6 K 
were compared to those of Shapiro and Moore (1976). The cooling function A 
was similar in the two cases. The major differences occurred .in the ionized 
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fractions at low temperatures (T < 3 x 10 K). These differences are probably 
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due to differences in the input atomic physics. In our caicuJation, the 

isobaric function A(T) was found to be close to the isochoric function. 

The ions of interest here have non-negj.iR.ible values of f ^ in the 

temperature range 10 4 - 10^ K. We found that a cooling gas with T > 1() R K 

initially remains in approximate ionization equilibrium until the temperature 
0 

drops to about 10 K. We also found that the time to reach ionization 

g 

equilibrium at T > 10 K for an initially underionized gas is shorter than the 

0 

cooling time of the hot gas. Our calculations therefore have T = 10 K and 

max 

4 

T . =10 K, and assume equilibrium ionization at T The initial density 

min max 

is not important for the ionization calculations because two body processes 

dominate both the ionization and recombination in the regime of interest so 

that the density scales out of the problem. The remaining parameter in 

equation (1) is (N^/n ). Column densities of selected ions are given in Table 
7 -1 

1 for N„/n = 10 cm s . The results can be scaled to other values for this 
H o 

quantity. The ions have been chosen because they have strong ultraviolet 

resonance lines (e.g. Savage 1984). They are C TV Ml 550. 8, 1548.2; 

Si IV 221402.8, 1393.8; N V 221242.8, 1238.8; S IV 21062.7; 0 VI 221037.6, 

1031.9; and S VI 22944.5, 933.4. The quantity T is the temperature at which 

a transition is made from isobaric (T > T ) to isochoric (T < T. ) evolution. 

tr tr 

The last line in Table 1 gives column densities that would be obtained using 

ionization equilibrium values for f and A, and assuming isobaric evolution. 

The column densities for higher ion stages (0 VI and N V among those 

quoted here) depend to some extent on the choice of T . in these cases f. 

max i 

0 

is small but not negligible at 10 K, so that changing T ^ to a value 
0 

> 3 x io K (where the f are negligible) can increase the column densities by 
factors of 2.2 (for 0 VI) and 1.6 (for N V). 
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The time-dependent case leads to the persistence of highly ionized 
species to lower temperatures than would be expected from equilibrium 
calculations. This is particularly true of Li-like species (e.g. C IV. N T V, 
and 0 VI) as the He-like ion stages recombine especially slowly. It can also 
be seen that the isobaric case gives lower column densities than the isochoric 
case. This is related to the property that n « 1/T during isobaric evolution, 
which increases the (n/n^) factor in equation (t). This effect is largest, for 
ions formed at low temperatures. It can be seen from Table 1 that the 
differences are larger for the lower ionization stages. 

The choices of the parameter T correspond to a selection of physical 
distance scales. If a region which has diameter D at T = 10 ' K is cooling 
from high temperatures, it will remain roughly isobaric until the cooling time 
becomes shorter than the sound crossing time. The transition temperature T 

5 

is the temperature at which this equality holds. For T = 4.7 x 10' K , we 

-1 5 

find that D ~ 2.8 n _ kpc , while for T x = 2.2 x 10 K, the appropriate value 

o-3 tr 

is D ~ 200 n * pc. where n 0 represents the initial H density in units of 
o-3 o-3 

_3 _3 

10 cm The distance scale D varies quickly with the transition 
temperature, [D « T^^ 6 /A(T)] so transition temperatures of a few xlo° K will 
be appropriate for many distance scales in an astrophysical context. For 
larger cooling regions, pure isochoric flow will be applicable, while for 
smaller ones, pure isobaric cooling will apply. 

In a similar way, we can calculate the strengths of emission lines that 
should be seen from these cooling plasmas. We find 


line 


= k 


M 

m. 


f n } 

i/' 

l HJ 


max 

/ 

i 

min 


(3 

12 


, f . j . . dT 
1 l line 

- 1 -v — — - ■ — 


( 2 ) 



7 


-3 - 1 

where j. . is the line emissivity (in erg cm s ), and L,. the emission 
lane line 

line luminosity (in units of erg s ^ ) . Table 2 lists luminosities from 
various strong ultraviolet lines which we have calculated for both the 
isochoric and isobaric cases. Since these in general vary by less than a 

factor of two, the intermediate cases are not shown in the table. The values 

-1 4 6 

are normalized to M = 1 Mo yr and T . and T are 10 K and 10 K, 

u min max 

respectively. 


III. APPLICATIONS 


i) The Galaxy 


Chevalier and Oegerle (1979) considered a model for the Galactic corona 

in which approximately 4 Mg yr 1 is circulated on each side of the galactic 

-3 -3 

plane and the typical H density is 10 ' cm . Another estimate of M can be 
obtained by assuming that supernova remnants evolve as adiabatic blast waves 
in a uniform medium. Then 


M = 23 f — 1 


- A f * H T 1 1 -.1 

(.40 yrj ^ 51 ^6 ’ M 0 Y r 

1 yrj ( 1 0° 1 erg J lio K J 


(3) 


is the rate at which mass is heated above temperature T , where z is the time 
between supernovae and E is the energy per supernova. If the supernova energy 


is eventually radiated, this gives an estimate of the mass cooling rate. 

-1 -3 -3 

If we take M = 4 Mg yr on one side of the plane, n Q = 30 cm and the 

2 ‘ 7-1 

area of the Galaxy to be xR with R - 15 kpc, then N ,/n ~ 1.6 x 10 cm s 

H o 

cj 5 

The relevant rows of Table 1 are the T. = 2.2 x 10' K and T. = 4.7 x 10' K 

tr tr 


ones, depending on the sizes of the cooling regions. The column densities in 
Table 1 must be multiplied by 1.6 to obtain the predicted column densities for 
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the GaJactic corona. The accuracy of this scaling factor is no better than a 
factor of 2; the accuracy within a row depends on uncertainties in the atomic 
physics and in the elemental abundances. 

3 3 -2 

The model column densities are thus N(C IV) = (4. 3-7. 9) x 10 ' cm 

N ( 0 VI) = (5. 8-6.0) x 10 14 cm” 2 . N(N V) = (2. 8-3. 6) x 10 13 cm” 2 , and N(S.i IV) 
12 -2 

= (3. 3-6. 4) x 10 cm Observed values perpendicular to the galactic plane 

are N(C IV) = 10 14 cm” 2 , N(N V) = 3 x io 13 cm” 2 , and N(Si IV) = 

13 -2 

2.5 x io cm (Pettini and West 1982; Savage and Massa 1985). The cooling 
model predicts value for N(N V) that is in excellent agreement with the 
observations. In contrast, photoionization model values for N(N V) typically 
fall short of the observed values by about an order of magnitude (e.g. 

Fransson and Chevalier 1985). On the other hand, photoi oni zation models 
predict the observed amount of Si IV while the cooling model underproduces 
Si IV by almost an order of magnitude. C IV can be produced in substantial 
quantity in either model, although the cooling model value may fall short of 
the observed value. 

As expected, the cooling model predicts a substantial column density of 

0 VI. There are no data on 0 VI in the halo, but it has been observed locally 

13 -2 

in the disk with a column density through the disk of about 2 x 10 cm 

(Jenkins 1978). This column density is considerably smaller than the model 

value obtained here; column densities of C IV and Si IV through the disk 

(Savage 1984) are also smaller than our model values. This may imply that 

supernova heated gas does not cool in the galactic disk (see also Cox 198.1). 

Another possibility is that the cooling occurs under high pressure conditions 

so that n for the cooling gas is high. For example, if supernova remnants 
o 

-8 

expand into gas with hydrogen density of about 0 . 2 cm , the shocked gas would 
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_3 

have a density close to 1 cm , and the column densities would be 
correspondingly decreased. 

Feldman, Brune, and Henry (1981) claim to have detected ultraviolet 

emission lines from the interstellar medium. They find fluxes of 9500 ± 3500 

(Si IV ,31397 + 0 IV] ,11406 ) , 10,000 ± 4000 (N IV ,11488), 8000 ± 3000 (C IV 

,11549), and 17,000 ± 4000 (0 III ,11663) in units of photons cm 2 s sr * . In 

the cooling model discussed above, C IV ,11549 (summing both members of the 

doublet) is the strongest line in the wavelength region of interest and it has 

-2 -1 -1 

a predicted flux of 890 photons cm s sr for isobaric evolution (see 
Table 2). This is an order of magnitude smaller than the claimed observation 
and, while there may be ways of increasing the prediction, it is necessary not 

to violate the limit on N(C IV). Another problem with the claimed detections 

is that N IV 31488 is observed to be comparable in strength to C IV ,11549, 
while it is predicted to be considerably weaker. Thus, while models for the 
ultraviolet emission lines can be found (e.g. Paresce et al . 1983), the most 
plausible cooling model is not consistent with the claimed detections. 

ii) Galaxies and Clusters of Galaxies 

Other spiral galaxies may have properties similar to our Galaxy, so that 
the above model approximately applies. Except for observations of supernovae, 
we have no information on column densities of highly ionized atoms in external 
galaxies. Absorption lines of C IV ,31549, Si IV ,31397 and possibly N V ,33241 
were observed in the ultraviolet spectrum of SN1979c in M100 (Panagia et al . 
1980), but it was not possible to clearly resolve the Galactic and M100 
features. Deharveng et al. (1986) searched for ultraviolet emission lines 
from a hot gaseous halo in the edge-on galaxy NGC 4244. An upper limit of 
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_0 _2 _ | 

about 10 erg cm s sr was set on any emission line. This is a factor 
of > 10 larger than the prediction for the strongest line in the standard 
cooling model. 

Some galaxies may have substantially smaller heavy element abundances 

than the 'cosmic' values used in § II. From eqn. (1), it can be seen that a 

decrease in n /n is compensated for by a decrease in A provided that emission 
L H 

from heavy elements dominates the cooling and that the relative abundances of 
the heavy elements remain constant. However, the increased cooling time 
allows the ionization to approach collisional equilibrium and the column 
densities are correspondingly decreased. 

Early-type galaxies and clusters of galaxies are thought to have cooling 
inflows that are roughly spherically symmetric. Mass flow rates are inferred 
to be 0.02 - 3 Mq yr 1 in elliptical galaxies (Thomas et a] . 1986) and up to 
400 Mq yr 1 in clusters of galaxies (Fabian, Nulsen. and Canizares 1984). The 
luminosities in various ultraviolet lines can be found directly from Table 2 
with the appropriate M scaling. Models of steady cooling flows (Mathews and 
Bregman 1978; White and Chevalier 1984) show that the gas cools through the 
temperature range of interest here over a small range in radius. Thus we find 


N H M 

— _ — - _ v 

n . u 2 an 

o 47tK p 

o 


( 4 ) 


where R is the radius at which the cooling occurs and v. is the infail 

m 

velocity at radius R. Models show that v. varies over a relatively small 

in 

range. For elliptical galaxies with core radii in the range 10-500 pc. White 


and Chevalier (1984) found v^ in the range 150 - 300 km s For cluster 
cooling flows with M of 3 to 100 Mq yr *, Mathews and Bregman (1978) found v. 


in 
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of 350 to 550 km s * . The column densities in Tabic 1 shouid be multiplied by 
a factor in the range 1.5 - 5.5 for these cooling flows. A significant 
difference from the Galactic case is that the densities in the temperature 
regime of interest can be large, e.g. n Q = 1 cm for T = 10 K in the 
cluster models of Mathews and Bregman (1978). Under these circumstances, the 
cooling is likely to be isocnoric in the region of interest. 

These results show that the detection of ultraviolet absorption lines in 
cooling flows is quite plausible if there is a suitable background continuum 
source. Cluster flows with nuclear continuum sources include Virgo, A1795, 
and Perseus (Hu, Cowie, and Wang 1985). Observations of the ions discussed 
here will cover a temperature regime intermediate between that covered by X- 
ray emission and that covered by optical emission. At present there is a 
discrepancy between the small cooling region indicated by optical observations 
and the large cooling regions indicated by X-ray observations (Hu, Cowie, and 
Wang 1985). Ultraviolet observations should help to clarify this discrepancy. 

This research was supported in part by NASA grant NAGW-764 and NSF grant 


AST-84 . 
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TABLE 1 


-2 7 

Column densities, cm , scaled to N„/n =1 . 0x10 cm s 

H O 








C IV 

N V 

0 VI 

Si IV 

S' IV 

S VI 


isochoric 

1.14(14) 

4.68(13) 

6.54(14) 

9.90(12) 

1.61(13) 

7.12(12) 

=4.7x 10 5 K 
tr 

4.94(13) 

2.24(13) 

3.63(14) 

4.04(12) 

6.29(12) 

4.01(12) 

T =2.2x10° K 
tr 

2.68(13) 

1.75(13) 

3.73(14) 

2.09(12) 

2.99(12) 

3.57(12) 

isobar ic 

1.36(13) 

1.64(13) 

3.83(14) 

4.14(11) 

1.36(12) 

3.29(12) 

equi librium 

1.47(13) 

1.54(13) 

3.69(14) 

3.38(11 ) 

1 .27(12) 

3.27(12) 
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TABLE 2 

UV Emission Line Luminosities: M = 1 Mq yr 

38 

(1.0 x 10 erg s ) 


species wavelength isochoric isobaric 

(A) 


Ne 

VI 

1006 

.0 

0, 

.214 

0 , 

.346 

Ar 

VI 

3008. 

,0 

0. 

193 

0 . 

317 

0 

VI 

1031 

.9 

7, 

.11 

11 . 

.0 

0 

VI 

1037 

.6 

3, 

.55 

5. 

.50 

S 

III 

1198 

.0 

0 

.300 

0 

.493 

S 

V 

1199 

.0 

0. 

,023 

0 . 

,039 

Si 

III 

1206 

.5 

0. 

.271 

0. 

.517 

0 

V 

1218. 

.0 

0. 

816 

1 . 

37 

N 

V 

1238 

.8 

0. 

.524 

0. 

.805 

N 

V 

1242 . 

8 

0. 

262 

0 . 

402 

Si 

IV 

1393 

.8 

0. 

.053 

0. 

.082 

0 

IV 

1402, 

.0 

0. 

,590 

0 . 

974 

Si 

[V 

1402 

.8 

0. 

.026 

0. 

.041 

C 

IV 

1548. 

.2 

0. 

955 

1 . 

52 

C 

IV 

1550 

.8 

0, 

.478 

0. 

.759 
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